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The microscopic approach to calculating the free energy of a three-dimensional 
Ising-like system in a homogeneous external field is developed in the higher non- 
Gaussian approximation (the p 6 model) at temperatures above the critical value 
of T c (T c is the phase-transition temperature in the absence of an external field). 
The free energy of the system is found by separating the contributions from the 
short- and long-wave spin-density oscillation modes taking into account both the 
temperature and field fluctuations of the order parameter. Our analytical calcula- 
tions do not involve power series in the scaling variable and are valid in the whole 
field-temperature plane near the critical point including the region in the vicinity 
of the limiting field h c , which divides external fields into the weak and strong ones 
(i.e., the crossover region). In this region, the temperature and field effects on the 
system are equivalent, the scaling variable is of the order of unity, and power series 
are not efficient. The obtained expression for the free energy contains the leading 
terms and terms determining the temperature and field confluent corrections. 

PACS numbers: 05.50. +q, 05.70.Ce, 64.60.Fr, 75.10.Hk 

1 Introduction 

Despite the great successes in the investigation of three-dimensional (3-D) 
Ising-like systems made by means of various methods (see, for example, [1]), 
the statistical description of the critical behavior of the mentioned systems 
in terms of the temperature and field variables and the calculation of scaling 
functions are still of interest [2] . The Ising model is widely used in the theory 
of phase transitions for the study of the properties of various magnetic and 
non- magnetic systems (ferroelectrics, ferromagnets, binary mixtures, etc.). 
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The description of the phase transitions in the 3D magnets, usually, is 
associated with the absence of exact solutions and with many approximate 
approaches for obtaining different system characteristics. In this article, the 
behavior of a 3D Ising-like system near the critical point in a homogeneous 
external field is studied using the collective variables (CV) method [3-5]. The 
main peculiarity of this method is the integration of short-wave spin-density 
oscillation modes, which is generally done without using perturbation theory. 
The CV method is similar to the Wilson non-perturbative renormalization- 
group (RG) approach (integration on fast modes and construction of an 
effective theory for slow modes) [6-8]. The term collective variables is a 
common name for a special class of variables that are specific for each indi- 
vidual physical system [3,4]. The CV set contains variables associated with 
order parameters. Because of this, the phase space of CV is most natural 
for describing a phase transition. For magnetic systems, the CV pk are the 
variables associated with modes of spin-moment density oscillations, while 
the order parameter is related to the variable p , in which the subscript "0" 
corresponds to the peak of the Fourier transform of the interaction potential. 

The free energy of a 3D Ising-like system in an external field at tem- 
peratures above T c is calculated using the non-Gaussian spin-density fluctu- 
ations, namely the sextic measure density. The latter is represented as an 
exponential function of the CV whose argument includes the powers with the 
corresponding coupling constants up to the sixth power of the variable (the 
p 6 model). 

The present publication supplements the earlier works [9-12], in which 
the p 6 model was used for calculating the free energy and other thermody- 
namic functions of the system in the absence of an external field. The p 6 
model provides a better quantitative description of the critical behavior of 
a 3D Ising-like magnet than the p 4 model [11]. For each of the p 2m models, 
there exists a preferred value of the RG parameter s = s* (s* = 3.5862 for 
the p 4 model, s* = 2.7349 for the p 6 model, s* = 2.6511 for the p 8 model, and 
s* = 2.6108 for the p 10 model) nullifying the average value of the coefficient 
in the term with the second power in the effective density of measure at the 
fixed point. The values of s close to s* are optimal for the given method of 
calculations. The difference form of the recurrence relations (RR) between 
the coefficients of effective non-Gaussian densities of measures operates suc- 
cessfully just in this region of s. It was established (see, for example, [11,12]) 
that as the form of the density of measure becomes more complicated, the 
dependence of the critical exponent of the correlation length v on the RG 
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parameter s becomes weaker gradually, and, starting from the sextic density 
of measure, the value of the exponent v, having a tendency to saturate with 
increasing m (which characterizes the order of the p 2m model, m = 2,3, 4, 5), 
changes insignificantly. The Ising model corresponds to the p 2m model ap- 
proximation, where the order of the model 2m > 4. The p 4 model allows us 
to go beyond the classical analysis and to describe all qualitative aspects of 
the second-order phase transition. The critical behavior of a 3D Ising-like 
system within the CV method can be described quantitatively at 2m > 6, 
and, in particular, at 2m = 6. It was shown in [11, 13] that the graphs of the 
temperature dependences of the order parameter (the spontaneous magne- 
tization) and specific heat for the p 6 model agree more closely with the Liu 
and Fisher's results [14] than the corresponding plots for the p 4 model. The 
correctness of the choice of the p 6 model for investigations is also confirmed 
in [15] and [16], where the effective potential is studied for the scalar field the- 
ory in three dimensions in the symmetric and spontaneously broken phases, 
respectively. In this case, probability distributions of average magnetization 
in the 3D Ising model in an external field, obtained with the help of the 
Monte Carlo method, were used. Tsypin [15, 16] proved that the term with 
the sixth power of the variable in the effective potential plays an important 
role. 

The methods existing at present make it possible to calculate universal 
quantities to a quite high degree of accuracy (see, for example, [1]). The 
advantage of the CV method lies in the possibility of obtaining and analysing 
thermodynamic characteristics as functions of the microscopic parameters of 
the initial system [11-13]. The results of calculations for a 3D Ising system 
on the basis of the p 4 and p 6 models are in accord with the results obtained 
by other authors (see [11,13]). In [17], the scaling functions of the order 
parameter and susceptibility, calculated on the basis of the free energy for the 
p 4 model, were graphically compared with other authors' data. Our results 
accord with the results obtained within the framework of the parametric 
representation of the equation of state [18] and Monte Carlo simulations [2]. 

The expressions for the thermodynamic characteristics of the system in 
the presence of an external field have already been obtained on the basis of 
the simplest non-Gaussian measure density (the p 4 model) in [19-22] using 
the point of exit of the system from the critical regime as a function of the 
temperature (the weak- field region) or of the field (the strong- field region). 
In [19,20], the thermodynamic characteristics are presented in the form of 
series expansions in the variables, which are combinations of the tempera- 



3 



ture and field. Our calculations in the p 4 model approximation were also 
performed for temperatures T > T c [21] and T < T c [22] without using sim- 
ilar expansions for the roots of cubic equations appearing in the theoretical 
analysis. In this article, the free energy of a 3D uniaxial magnet within the 
framework of the more complicated p 6 model is found without using series 
expansions introducing the generalized point of exit of the system from the 
critical regime. This point takes into account the temperature and field vari- 
ables simultaneously. In our earlier article [17], the point of exit of the system 
from the critical regime was found in the simpler non-Gaussian approxima- 
tion (the p 4 model) using the numerical calculations. In contrast to [17], the 
point of exit of the system in the present article is explicitely defined as a 
function of the temperature and field. This allows one to obtain the free en- 
ergy without involving the numerical calculations that is our problem solved 
in the present article. 



2 Integration of partition function of the sys- 
tem in the p 6 model approximation taking 
into account effect of an external field 

We consider a 3D Ising-like system on a simple cubic lattice with N sites 
and period c in a homogeneous external field h. The Hamiltonian of such a 
system has the form 

H = -\ E^jO^i-^E^ (i) 

where is the distance between particles at sites j and 1, and <7j is the 
operator of the z component of spin at the jth site, having two eigenvalues 
+ 1 and -1. The interaction potential is an exponentially decreasing function 

$(r jI )=Aexpf-?V (2) 



b 

Here A is a constant and b is the radius of effective interaction. For the 
Fourier transform of the interaction potential, we use the following approxi- 
mation [3, 11, 12]: 

nk) - jo, B'<k<B, {6} 
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where B is the boundary of the Brillouin half-zone (B = tt/c), B' = (by/2)- 1 , 
<£(0) = 8nA(b/c) 3 . 

In the CV representation for the partition function of the system, we 
have [3,23] 



exp 



J(p) (dpf. 



(4) 



Here the summation over the wave vectors k is carried out within the first 
Brillouin zone, (3 = l/(kT) is the inverse thermodynamic temperature, the 
CV pk are introduced by means of the functional representation for operators 
of spin-density oscillation modes pk = (y/N)' 1 ci ex P( — zkl), 



J(P) 



2 N I exp 



2m E^kPk + E( 2 ™) N ' 

k n>l 



M 2n \ - , 

\ Ln )- k lv ..,k 2 „ 



(dbJ 



(5) 



is the Jacobian of transition from the set of N spin variables <7\ to the set 
of CV pk, and <5k 1 +...+k 2 „ * s the Kronecker symbol. The variables Wk are 
conjugate to pk, and the cumulants M.2n assume constant values (see [3-5]). 

Proceeding from Eqs. (TJ]) and ([5]), we obtain the following initial expres- 
sion for the partition function of the system in the p 6 model approximation: 



Z 



2 N 2 (N'-l)/2 e a' N' 



exp 



-a[(N') 1/2 Po 



l 2l 



k 

k<B' 



1=2 



X 



E Pk x • ■ ■ Pk 2! 5ki+-+k 2 



kl,...,k 2i 



(dp) 



N' 



(6) 

7iV2b/c, 



Here iV = iVso d (d = 3 is the space dimension), s = B/B' 
and a ; = -s d Q /2 h', ti = (3h. The expressions for the remaining coefficients 
are given in [9-12]. These coefficients are functions of so, i.e., of the ratio of 
microscopic parameters b and c. The integration over the zeroth, first, second, 
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. . ., nth layers of the CV phase space [3-5, 11] leads to the representation 
of the partition function in the form of a product of the partial partition 
functions Q n of individual layers and the integral of the "smoothed" effective 
measure density 



2 N 2^- i y 2 Q Q 1 --.Q n [Q(P n )] N ^ 



(7) 



The expressions for Q n , Q(P n ) are presented in [9-12], and N n+ i = 
N f s -d(n+i)_ The gextic 

measure density of the (n + l)th block structure 

M n+1 \p) has the form 



(n+l) 



(P) 



= exp 



(n+l) N l/2 
- a l iV ra+lP0 



2 z ^ 



d n +i(k)pkP-k 



k 

k<B n+1 



a 



(n+l) 



21 



3 

1=2 (2l) } -N l n+ \ kl ,...,k 2i 

k i< B n+l 



J2 fa 



■pk 2l 



•+k 2i 



(8) 



(n+l) 



a 2 ' — f3$(k), di +V> and a^ +1 ^ are the 



where B n+1 = B's^ n+1 \ d n+1 (k) 
renormalized values of the coefficients a[ and a' 2Z after integration over n+l 



layers of the phase space of CV. The coefficients = s n t n , d n (0) = s 



-2n„ 



» 



~ 4n u n , and 



s 6n w r 



[appearing in d n (k) = d n (0) + 2(3®(0)b 2 k 2 }, a!£ 
are connected with the coefficients of the (n + l)th layer through the RR 



r n +\ = s 



Un+l ■ 
Wn+l 



q + uj Y(h n , a r 

u n B(h n ,a n ), 
s^ d u^D{h n ,a n ) 



A-d 



(9) 



whose solutions 



TV, 



= r 



(0) +CljE n + C2W (0) (M (0) r l/ 2jE n 



+c 3 ^gV {0) 



1 T7<n 
4 ) 
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Un = u^+c lW { §{u^ 2 E^ + c 2 E n , 

Wn = wW + cxw^u^EZ 

+c 2 wg\u^ 2 E^ + c 3 E2 (10) 

in the region of the critical regime are used for calculating the free energy of 
the system. Here 

Y(h n ,a n ) = s d/2 F 2 ( Vn ^ n ) [C(h n ,a n )}' 1/2 , 

B(h n , a n ) = s 2d C(i] n , £ n ) [C(hn, a„)] _1 , 

D(h n , a n ) = s 7d ^N( Vn , U) [C(K, a n )}- 3/2 . (11) 

The quantity q = q/3$>(0) determines the average value of the Fourier trans- 
form of the potential /3<&(B n+ i, B n ) = /3$(0) — q/s 2n in the nth layer (in this 
article, q — (1 + s~ 2 )/2 corresponds to the arithmetic mean value of k 2 on 
the interval (1/s, 1]). The basic arguments h n and a n are determined by the 
coefficients of the sextic measure density of the nth block structure. The 
intermediate variables r] n and £ n are functions of h n and a n . The expressions 
for both basic and intermediate arguments as well as the special functions 
appearing in Eqs. (jTTI) are the same as in the absence of an external field 
(see [9-12]). The quantities E[ in Eqs. (fTDl) are the eigenvalues of the matrix 
of the RG linear transformation 





t n +l 


_ t (0) 




( 


R n 













( 


tfl 


_f(0) 






r n +i 


_ r (0) 









R-22 


-R23 


R24 




r n 


_ r (o) 




U n +1 


_ u (0) 









R32 


-R33 


R34 








_ u (0) 




V 


W n +1 


- w (0) 


) 


V 





R42 


-R43 


R44 


) 


\ 


W n 


- w (f>) 


) 



(12) 



We have E\ = Rn = s^ d+2 ^ 2 . Other nonzero matrix elements Rij (i = 2, 3, 4; 
j = 2,3,4) and the eigenvalues E 2 , E 3 , E4 coincide, respectively, with the 
quantities R^^ (i\ = i — 1; ji = j — 1) and E ly E 2 , E 3 obtained in the 
case of h = 0. The quantities /o, ipo, and i[)q characterizing the fixed-point 
coordinates 

t (o) = 0> r (o) = _ /o/3 <£( ), 

«w = MPm) 2 , = MPm? (13) 

as well as the remaining coefficients in Eqs. ( fl~0l) are also defined on the basis 
of expressions corresponding to a zero external field. 
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3 Using the generalized point of exit of the 
system from the critical-regime region for 
calculating the free energy 

Let us calculate the free energy F = —kTlnZ of a 3D Ising-like system 
above the critical temperature T c . The basic idea of such a calculation on 
the microscopic level consists in the separate inclusion of the contributions 
from short-wave (Fqr, the region of the critical regime) and long-wave (Flgr, 
the region of the limiting Gaussian regime) modes of spin-moment density 
oscillations [3-5]: 

F = F + F CR + F LGR . (14) 

Here Fq = —kTN In 2 is the free energy of iV noninteracting spins. Each 
of three components in Eq. (}Tlj) corresponds to individual factor in the 
convenient representation 

Z = 2 Z C rZ LG r (15) 

for the partition function given by Eq. (J7]). The contributions from short- 
and long-wave modes to the free energy of the system in the presence of 
an external field are calculated in the p 6 model approximation according to 
the scheme proposed in [9-12]. Short-wave modes are characterized by a 
RG symmetry and are described by the non- Gaussian measure density. The 
calculation of the contribution from long-wave modes is based on using the 
Gaussian measure density as the basis one. Here, we have developed a direct 
method of calculations with the results obtained by taking into account the 
short-wave modes as initial parameters. The main results obtained in the 
course of deriving the complete expression for the free energy of the system 
are presented below. 

3.1 Region of the critical regime 

A calculation technique based on the p 6 model for the contribution Fqr 
is similar to that elaborated in the absence of an external field (see, for 
example, [5,10,11]). Carrying out the summation of partial tree energies F n 
over the layers of the phase space of CV, we can calculate F C r. 

FcR = K + F CR> 
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Fq = -kT N' [In Q(M)+\nQ(d)l 

rip 

F' C R= E F n- (16) 



n=l 



An explicit dependence of F n on the layer number n is obtained using solu- 
tions fflUj) of RR and series expansions of special functions in small deviations 
of the basic arguments from their values at the fixed point. The main pecu- 
liar feature of the present calculations lies in using the generalized point of 
exit of the system from the critical regime of order-parameter fluctuations. 
The inclusion of the more complicated expression for the exit point (as a 
function of both the temperature and field variables) [24] 

= _ Mh' + hp _ 

' 21n£, [ ' 

leads to the distinction between formula ([TBI) for F' CR and the analogous 
relation at h = [10,11]. The quantity h = h'/f is determined by the 
dimensionless field h f , while the quantity h c = f po is a function of the reduced 
temperature r = (T — T c )/T c . Here f = cf r/fo, po = ln^i/lni^ = 
(d + 2)u/2, Ci characterizes the coefficient c\ in solutions fflUj) of RR, v = 
In s/lnE2 is the critical exponent of the correlation length. At h = 0, n p 
becomes m T = — \nf/\nE 2 — 1 (see [5,10,11]). At T = T c (r = 0), the 
quantity n p coincides with the exit point = — In h/ hiEi — 1 [25]. The 
limiting value of the field h c is obtained by the equality of the exit points 
defined by the temperature and by the field (m T = n^). 

Having expression (fT7|l for n p , we arrive at the relations [26] 



E^ +1 = {h 2 + hl)- l/ \ fE^ +1 = H c , 
H e = hl'*(h 2 + hl)- 1 K 2po \ 
E n z p+1 = H 3 , H 3 = {h 2 + ~h 2 c ) A ^°\ 
E2 P+1 = Hi, H 4 = (h 2 + ~h 2 c ) A ^ 2po \ 

s -(n p+ l) = { J l 2 + ~ h 2 ) l/(d+2)^ (lg) 

where A\ = — InE^/ lni?2 and A2 = — In £4/ In E2 are the exponents, which 
determine the first and second confluent corrections, respectively. Numerical 
values of the quantities Ei {I = 1, 2, 3, 4), z/, A 1; and A 2 for the optimal RG 
parameter s = s* = 2.7349 are given in Table [H In the weak-field region 
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Table 1: The eigenvalues Ei and the exponents is, Ai, A2 for the p 6 model. 



Ei E 2 


E 3 


£4 


v Ai A 2 


12.3695 4.8468 


0.4367 


0.0032 


0.637 0.525 3.647 



(/i <C /i c ), quantities (ITS1) can be calculated with the help of the following 
expansions: 

E^ +l = h~ 1 (l-- 1 t), h- 1 = r- po , 
1 c \ 2 h?J' 

1 h 2 



H c =l 

2p h 2 r 



c 

"2" 



-(np+l) = / j + 1 £\ 

In the strong-field region (h 3> h c ), these quantities satisfy the expressions 



f 3 = f 1 + ^ 



s -(n P+ l) = ~ h 2/ {d+ 2) ( ^ _J_ K\ (2Q) 

V d + 2 h? ) K ' 
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It should be noted that the variables h/h c (the weak fields) and (h c /h) l l po 
(the strong fields) coincide with the accepted choice of the arguments for 
scaling functions in accordance with the scaling theory. In the particular case 
of h = and r ^ 0, Eqs. (pi| are defined as £™ p+1 = f- p \ H c = l,H 3 = f A \ 
H 4 = f A2 , s~^ +1 ^ = r v . At h ^ and r = 0, we have ~hE™ p+1 = 1, H c = 0, 
H 3 = h A ^ Po , H 4 = h A */P°, = ^2/(^+2) [ see Eqs _ p}|)]. 

We shall perform the further calculations on the basis of Eqs. ffTSj) . 
which are valid in the general case for the regions of small, intermediate 
(the crossover region), and large field values. The inclusion of E 3 P+1 (or H 3 ) 
leads to the formation of the first confluent corrections in the expressions for 
thermodynamic characteristics of the system. The quantity E"2 P+1 (or H4) 
is responsible for the emergence of the second confluent corrections. The 
cases of the weak or strong fields can be obtained from general expressions 
by using Eqs. ffT9l) or ([20]) . We disregard the second confluent correction 
in our calculations. This is due to the fact that the contribution from the 
first confluent correction to thermodynamic functions near the critical point 
(r = 0, h = 0) for various values of s is more significant than the small 
contribution from the second correction (h 2 + h 2 . -C 1, is of the order of 
0.5, and A 2 > 2, see Tableland [9]). 

Proceeding from an explicit dependence of F n on the layer number n 
[5,9,10] and taking into account Eqs. (Fl8|) . we can now write the final 
expression for F C r ffTEj) : 



F CR = -kTN' ( 7 f R) + 7l r + 72 r 2 ) + F s , 



(21) 



Here characterizes c 2 in solutions (FIU|) of RR, 

,(0) ,(1) -1/2, JT 

-(CR)(0)+ JCR , JCRfo JOMc 

73 - 1 — a + 



-(CR)(1) 



73 = 1 ET— 5 + 



1-s" 3 1-E 2 s- 3 
1 - Efs" 3 ' 

,(2) 1 ,(4) -3/2, rr 

JCrVo . JcRfo J0 H c 



1-E 3 s- 3 l-E 2 E 3 s- 3 

1 fcR ( Po 2 (foHc) 2 / 99 s 
f l-E 2 E 3 s-z ' 1 j 
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and the coefficients 



Ik = iSP + $\ fc = l,2 (23) 
are determined by the components of the quantities 

(0) , (1) , (2) 2 

7o = 7o +7o T + To T , 

* = + tfV + *SV. (24) 

The components 5q (« = 0, 1, 2) satisfy the earlier relations [5,9,10] obtained 
in the case of a zero external field. The components 7g are given by the 
corresponding expressions at h — under condition that the eigenvalues E%, 
E 2 , and E 3 should be replaced by E 2 , E 3 , and E4, respectively. 

Let us now calculate the contribution to the free energy of the system from 
the layers of the CV phase space beyond the point of exit from the critical- 
regime region. The calculations are performed according to the scheme pro- 
posed in [3,5,11,12]. As in the previous study, while calculating the partition 
function component Z LGR from Eq. (fT5l) . it is convenient to single out two 
regions of values of wave vectors. The first is the transition region (Z LGR ) 
corresponding to values of k close to B n , while the second is the Gaussian 

region (Z^ GR ) corresponding to small values of wave vector (k — > 0). Thus, 
we have 

Zlgr = Z L x GR Zl GR . (25) 
3.2 Transition region 

This region corresponds to m layers of the phase space of CV. The lower 
boundary of the transition region is determined by the point of exit of the 
system from the critical- regime region in = n p + 1). The upper boundary 
corresponds to the layer n p + rh + 1. We use for m the integer closest to 
m' Q . The condition for obtaining m' is the equality [10, 11] 

I K p +m' |= : A °, , (26) 



where A is a large number (A > 10). 
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The free energy contribution 

IlgrA™) = ln (^) + \ In 24 - ^ \nC(r] np+m ,£ np+m ) 

+ ln Io{h np + m +l, «n p +m+l) 

(27) 

corresponding to Z^ 1 ^ from Eq. (I25I) is calculated by using the solutions of 
RR. 

The basic arguments in the (n p + m)th layer 

hn p +m = ( r n p +m + <?) (6/u„ p+m ) 1/<2 , 

^6 / 3/2 , 98 n 

ttn p + m — w n p +m/U np+m [ZO ) 

can be presented using the relations 

= -sffoEr'Hh 2 + hir x '\ 

r np+m = /3$(0)(-/ + /o^ 2 m - 1 

+C 20 -H 3 <£ W 12 ii 3 

S +m = (/5<I(0)) 2 U + foH^w^E™- 1 



I „(°) o" pm-1 



= (/?<l(0)) 3 (^o + /o^o^i ) ^ 2 m - 1 

+4?^3^ /2 ^° 2 ) ^- 1 ) (29) 



obtained on the basis of Eqs. (jTUl) and ffT8l) . We arrive at the following 
expressions: 



_ 1,(0) fi , IW JO) tr 
n n p +m — lb n p +m \ L > u n p +m c 20 n Z 



, (o) 9 — /o + foH c E™ 1 

i^n-p+m — V O 



/ -1/2 (0) 

1(1) _ pm-l( V?0 ™12 

' ' r) „ _|_<m -C^' 



\Q - Jo + Jo^cE 2 
13 



2 VQ + f Q H cV y 2 w^E^)' 



(o) /, 



a 



(0) 



+771^20 



(0) pm-1 



V6 t/> + f H c <p w$>E% 
15 (ipo + foH^w^ET 1 )^'' 



a 



(i) 



-^r 1 1 



1/2 (0) 
¥>0 ^32 



2 ^ + foH^w^ET 1 



(30) 



In contrast to if c , the quantity H 3 in expressions (1301) for h np+m and a np+m 
as well as in expression (I2T!) for F s takes on small values with the variation 
of the field h (see Fig. [T]). The quantity H c at h — > and near /i c is close to 




h/hc 



Figure 1: Dependence of quantities H c and if 3 on the ratio h/h c for the RG 
parameter s = s* = 2.7349 and the reduced temperature r = 10~ 4 . 

unity and series expansions in H c are not effective here. 

Power series in small deviations (h np + m — hnp+m) &nd (ct np + m — ^np+m) 
for the special functions appearing in the expressions for the intermediate 
arguments 

Vn p+m = (Qs d ) l ' 2 F 2 (h np _ 



+771 



X 

15" 



d/2 N(h np 



-1/2 



X 



\1 ~ 3 / 2 



(31) 
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allow us to find the relations 

Vrip+m 



„(0) 
'Irip+m 



1 - (t 

-(np+m) (0) (1) 



(n p +m),(0) 7(1) 

' ''np +m ""np+m 



4~np 



+m 



//■V " ""' 0>n p +mUn p +m) 4o #3 
1 



t(0) 
Srip+m 



cK+m)i.(0) 1(1) 
SI '''np+m'''np+m 



(n p +m) (0) 



ft. 



n p +m a n p +rn ) c 20 -"3 



(1) 



Con -"a 



(32) 



mi (0) _(n p +m) _(n p +m) , *(0) /.(n p +m) *(n„+m) t 

The quantities ?7 n ; +m , ?7i , rfo and G; +m , Ci , G are func- 



tions Of F 



*(n p +m) 



2/ 



r*(n p +m) i j*(n v +m) 
L 2l 



where 



j*(n p +m) 
L 2l 



'JO 

^ x exp( h np _^ m x x o; n?) _|_ m x ) dx. 



(33) 



Proceeding from expression (J2"T|) for fLGR^i™), we can now write the 
following relation accurate to within if 3 : 



(0) 

LGRi 



m 



LGRi (P^) 



(n„+m) 

(n p +m+l) 
^3 

(n p +m+l) 
(fi4 



111 



7T 



+ - In 24-^ InC^+^gJ 



.(0) 



+mj 



+ lnlo(fc 



(0) 

rip+m+l) u n p +m4- 



(0) 



-l) 



) 4"n p +m) 

(n p +m),(0) 7(1) 
Sri ' l 'n p +m''"n p +m 



(0) 



(np+m) (0) 



+^2 

(np+m+1), (0) 
-|-<^3 Ai 



ft 1 

n p +m+ 1 "tip +m+ 1 



(np+m+ l) (0) (i) 
-1-^4 

, (np+m) 



^np+m+1 ^n p +m+l J 



p(n p +m) 



/4, 



fc = l,2, 



-F 



-i*(n p +m+l) 
2 > 
*(n p +m+l) 



(34) 



The quantities &[™ p+m ' > , P 4 ^ p+m ' ) depend on p*( n p +m } 

J**K+'")// o **(n P +m) ) where 



as well as on F 



**(n p +m) 



2i 



T**(n p +m) 
l 2l 



J X exp(— Tjnp +m X — X — s^n.p+m^ 6 ) dx. 



(35) 
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The final result for F^q r [see Eqs. ( 12 7 jl and (134|) ] assumes the form 



LGR 



m 

JTR 



f(l) 



"to 



-3m AO) 
JLGRi 

-3m 7(1) 



3 > 



m=0 



m=0 



7. 



(36) 



On the basis of Eqs. (|26|) and (I30I) . it is possible to obtain the quantity m' 
determining the summation limit rfio in formulas 



m'r 



In L — In H c 



ln£ 2 +1 ' 

L = A 1 + (^-A 2 ) 1 / 2 , 



A 2 



!-■?■ + 



^WO w 21 



/o 12/ (1- S - 3 ) 2 



6/o 2 (l 



,-3^2 



(37) 



Let us now calculate the contribution to the free energy of the system 
from long-wave modes in the range of wave vectors 



k < B's- n 'p. 



n p + m + 2 



(38) 



using the Gaussian measure density. 



3.3 Region of small values of wave vector (k 

The free energy component 



?(2) 
LGR 



~-kT 



B n , 

X, VaP^ + ^lnd^k)-^ 



k=0 



dndO) 



o) 



(39) 



corresponding to Z L q R from Eq. (1251) is similar to that presented in [5,10,11]. 
The calculations of the first and second terms in Eq. (139]) are associated with 
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the calculations of the quantities 



2h n ' iF 2 (h n > i, a n / i) 



x 



d n '-i(B n ' , B n i i) 



-i 



P. 



K-i) 



-i 



where 



d n > p -i(B n , p , B^) = s ^'^{r^ + q), 

aadr^.i, h n<v _ x = h^_ ± (l + h^cffl H 3 ) , a^_ a = (l + 
satisfy the corresponding expressions from Eqs. (1291) and fl30|) at m = 
Introducing the designation 

P = h n > p -iF 2 {h n > p -i,a nlp _i) 

and presenting it in the form 

V 

we obtain the following relations for the coefficients: 



Po 

Pi 



a n'-lP20 



-1 



(n£-l) (0) _(1) 

+P22 a^-ia^-i- 



The quantities 



(n' p -l) _ *(n> p -l) (n' p -l) _ F 4 

P20 — r 2 1 P21 



*K-i) 



*K-i) 



K-i) 
P22 



8 



determine the function 



K-i) 
P20 



1 IP21 'V„-l n 'n'-l 



.^K-iUo) -(1) \ (0)rr 
+F22 a n'~l a n'-l C 20 U - 
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Here - r*^- 1 ) / T < n 'p- X ) 



i;r where 



r*K-l) 



L 2l 



x 21 exp(— h^/^x 2 — x 4 — a^_ 1 x ti ) dx. 



,(0) „2 _ „4 _ _(°) ^ 
i/ - 



(47) 



Taking into account Eqs. (14TI) and (|43|) . we rewrite formulas (I4TJ]) as 

P 2 K_1) = {i s - 2( ^- 1) /3<l(o)po(g-/o + /o^ ) 



1+ 



t-Pl I 4? if; 



v g-/o + /o^ 2 mo 

J n ^(A;) = S - 2{n ''>- 1) f3$(0)G + 2f3$(0)b 2 k 2 , 
G = g (l + g 1 c i SH 3 ), 

90 = \ [(-f + f H c E^)p + (p -2)q 

91 = ^ -\pi(q - fo + foH e E*°) 

1 9o 

The second term in Eq. ( |39i) is defined by the expression 



(48) 



52 In dn> p (k) = N K { \ ln(G + s" 2 ) + In s - nj, In s 
z fc=0 ^ 

+i ln(/?$(0)) - i + Gs 2 



-(G S 2 ) 3 / 2 arctan [(Gs 2 )" 1 / 2 ] }. 



(49) 



Relations (148 p and (I49p make it possible to find the component F^q r in 
the form 



p( 2 ) —JcT 
LGR ~ K1 



f(0)' _ s -3(m +l) jr(O) J(l)' = s -3(m +l) HI) 
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./ 



2 + 9o 



•<°> = -- In S - 

2 \ g + q 



-9o 



1 — J g' arctan 



9o9i 



2 U + ? (^o)- 1 + l 




-9o9i 



1 — - -y/gQ arctan 



= s 



2mo 



/(2flb). 



(50) 



#o = s 2 g , 74 

On the basis of Eqs. (155]) and (I50p . we can write the following expression 
for the general contribution F LGR = F^ GR + Fj? GR to the free energy of the 
system from long-wave modes of spin-moment density oscillations: 



H LGR 



= -kT 



mmi s2K+1)(1 _ ^ 

/J$(0) v 20 

Jlgr = frn + > / = 0, 1. 



(51) 



4 Total free energy of the system at T > T c 
as function of temperature, field and mi- 
croscopic parameters 

The total free energy of the system is calculated taking into account Eqs. 
(TH|) . (1211) . and (1511) . Collecting the contributions to the free energy from all 
regimes of fluctuations at T > T c in the presence of an external field and 
using the relation for s~( n p +1 ) from Eqs. ffl8l) . we obtain 



-kTN 



l'o + l[r + ^ + ^r + 4 1)+ c^H 3 



Ah' + ^) 3/5 + ^^(l-^4 0^3)(^ + hi)- 2 " 

7 ; = ln2 + So 3 7 r ) , 7i 



s o 3 7i, 7 2 = s 3 72, 
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# + = So- 3 (-7f H,( ' ,+ +m B ), 1 = 0,1. 



(52) 



The coefficients 7q , 7i, 72 are defined by Eqs. f|23l) . ^ is presented in 
Eqs. (HHI) . and 74" is given in Eqs. (!50|) . The coefficients of the non-analytic 
component of the free energy F [see Eqs. (1521) ] depend on H c . The terms 
proportional to H3 determine the confluent corrections by the temperature 
and field. As is seen from the expression for F, the free energy of the system 
at h = and f = 0, in addition to terms proportional to 

f 3u ( Qr ^6/5) and 

h 6 / 5 , contains the terms proportional to f 3u+A ^ and h 6 ^ 5+Al ^ Po , respectively. 
At h 7^ and f ^ 0, the terms of both types are present. It should be noted 
that Ai > Ax/pq. At h = h c , we have f Al = h Al ^ Po and the contributions 
to the thermodynamic characteristics of the system from both types of the 
corrections become of the same order. 

The advantage of the method presented in this article is the possibility of 
deriving analytic expressions for the free-energy coefficients as functions of 
the microscopic parameters of the system (the lattice constant c and param- 
eters of the interaction potential, i.e., the effective radius b of the potential, 
the Fourier transform $(0) of the potential for k — 0). 

5 Conclusions 

An analytic method for calculating the total free energy of a 3D Ising-like 
system (a 3D uniaxial magnet) near the critical point is developed on the 
microscopic level in the higher non-Gaussian approximation based on the sex- 
tic distribution for modes of spin-moment density oscillations (the p 6 model). 
The simultaneous effect of the temperature and field on the behavior of the 
system is taken into account. An external field is introduced in the Hamil- 
tonian of the system from the outset. In contrast to previous studies on the 
basis of the asymmetric p 4 model [19,20,27], the field in the initial process 
of calculating the partition function of the system is not included in the Ja- 
cobian of transition from the set of spin variables to the set of CV. Such 
an approach leads to the appearance of the first, second, fourth, and sixth 
powers of CV in the expression for the partition function and allows us to 
simplify the mathematical description because the odd part is represented 
only by the linear term. 

The theory is being built ab initio beginning from the Hamiltonian of the 
system up to the expression for the free energy. The main distinctive feature 
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of the proposed method is the separate inclusion of the contributions to the 
free energy from the short- and long-wave spin-density oscillation modes. 
The generalized point of exit of the system from the critical regime contains 
both the temperature and field variables. The form of the temperature and 
field dependences for the free energy of the system is determined by solutions 
of RR near the fixed point. The expression for the free energy F obtained 
at temperatures T > T c without using power series in the scaling variable 
and without any adjustable parameters can be employed in the field region 
near h c (the crossover region). The limiting field h c satisfies the condition 
of the equality of sizes of the critical-regime region by the temperature and 
field (the effect of the temperature and field on the system in the vicinity 
of the critical point is equivalent) [19,20,25,27]. In the vicinity of h c , the 
scaling variable is of the order of unity and power series in this variable are 
not effective. We hope that the proposed method as well as our explicit 
representations may provide useful benchmarks in studying the effect of an 
external magnetic field on the critical behavior of 3D Ising-like systems within 
the framework of the higher non-Gaussian approximation (the p 6 model). 
Proceeding from the expression for the free energy, which involves the leading 
terms and terms determining the temperature and field confluent corrections, 
we can find other thermodynamic characteristics (the average spin moment, 
susceptibility, entropy, and specific heat) by direct differentiation of F with 
respect to field or temperature. 
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